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ABSTRACT 

The probability distribution functions (PDF) of density of the interstellar medium (ISM) in galactic 
disks and global star formation rate are discussed. Three-dimensional hydrodynamic simulations show 
that the PDFs in globally stable, inhomogeneous ISM in galactic disks are well fitted by a single log- 
normal function over a wide density range. The dispersion of the log-normal PDF (LN-PDF) is larger 
for more gas-rich systems, and whereas the characteristic density of LN-PDF, for which the volume 
fraction becomes the maximum, does not significantly depend on the initial conditions. Supposing 
the galactic ISM is characterized by the LN-PDF, we give a global star formation rate (SFR) as a 
function of average gas density, a critical local density for star formation, and star formation efficiency. 
Although the present model is more appropriate for massive and geometrically thin disks (~ 10 pc) in 
inner galactic regions (< a few kpc), we can make a comparison between our model and observations 
in terms of SFR, provided that the log-normal nature of the density field is also the case in the real 
galactic disk with a large scale height (~ 100 pc). We find that the observed SFR is well-fitted by 
the theoretical SFR in a wide range of the global gas density (10 — 10 4 M Q pc -2 ). Star formation 
efficiency (SFE) for high density gas (n > 10 3 cm -3 ) is SFE = 0.001 — 0.01 for normal spiral galaxies, 
and SFE = 0.01 — 0.1 for starburst galaxies. The LN-PDF and SFR proposed here could be applicable 
for modeling star formation on a kpc-scale in galaxies or numerical simulations of galaxy formation, 
in which the numerical resolution is not fine enough to describe the local star formation. 

Subject headings: ISM: structure, kinematics and dynamics — galaxies: starburst — method: numer- 
ical 



1. INTRODUCTION 

The interstellar medium (ISM) in galaxies is characterized by highly inhomogeneous structure with a wide variety 
of physical and chemical states (|Mverslfi978| ). Stars are formed in this complexity through gravitational instability 
in molecular cores, but the entire multi-phase structures on a global scale are quasi-stable. It is therefore important 
to understand theoretically the structure of the ISM in a wide dynamic range to model star formation in galaxies. 
Observations suggest that there is a positive correlation between global star formation rate (SFR) and the average gas 
density: oc with N ~ 1.4 in nearby galaxies (|Kennicu t3[l99l) 3 . Since star formation process itself is a local 
phenomenon on a sub-parsec scale, the observed correlation between the structure of the ISM on a local scale and the 
global quantities, such as the average gas density, implies that the ISM on different scales is phy sically related. 

In fact, two- and three-dimensional hydrodynamic and magnetq - hydrodynamic s imulations fe.g. lBania fc Lvonlll98Cil : 
IVazquez-Semadeni. Gazol. fc Sca lo 2000: [Rosen fc Bregma n 1995; d e Avillezll2000f ) show that there is a robust relation 
between the local and global structures of the multi-phase ISM, which is described by a log-normal density PDF 
(probability distribution function). Elmegreen (2002) first noticed that if the density PDF is log- normal and star 
formation occurs in dense gases above a critical density, the Sch midt-Kennicutt law is repro duced. This provides 
a new insight on the origin of the scaling relation. More recently iKrumholz fc McKed (|2005h give a similar model 
on SFR in molecular clouds. In these theoretical predictions, dispersion of the LN-PDF, a, is a key parameter. 
Elme gr een (2002) used a = 2.4, w hich is taken from 2-D hydrodynamic simulations of the ISM (|Wada fc Norman] 

2001) . IKrumholz fc McKec (2005) assumed an empirical relation between a and the rms Mach number, which is 
suggested in numerical simulations of isothermal turbulence (see also §4.1). Therefore, it is essential to know whether 
the LN-PDF in galactic dicks is universal, and what determines a. However, most previous simulations, in which LN- 
PDF or power-law PDF are reported, are 'local' simulations: a patch of the galactic disk is simulated with a periodic 
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3 Recent observat i ons r eported a wide variety on the slope, e.g. N ~ 1.1 or 1.7 depending on the extinction models (Wong Sz Blitz 

2002) . Komugi et al. (2005) found that N ~ 1.33 for the central part of normal galaxies, and they also suggest that SFR is systematically 
lower than those in starburst galaxies (See also §3.3). 
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boundary condition. Apparently, such local simulations are not suitable to discuss statistical nature of the ISM in 
galactic disks. For ex ample, number of density condensations are not large enough (typically a few) (see IScalo et al.l 
Il998t ISIvz et~alll2005h . 

On the other hand, global hydrodynamic simulations for two-dimensional galactic disks or three-dimensional circum- 
nuclear gas disks suggested that the density PDF, especially a high-density part is well fitted by a single log-normal 
function over 4-5 decades (Wada & Norman 1999, 2001, hereafter WN01 and Wada 2001). The log-norm al PDF is 
also s een in a high-z galaxy formed by a cosmological N-body/AMR( Adaptive Mesh Refinement) simulation (|Kravtsovl 
2003). Nevertheless universality of LN-PDF and how it is related to global quantities are still unclear. 

In this paper, we verify the log-normal nature of the ISM in galactic disks, using three-dimensional, global hydrody- 
namic simulations. This is an extension of our previous two -dimensional studies of the ISM in galactic disks (WN01) 
and three-dimensional model in the galactic central region (|Wada fc Normanll2002l ). We confirm that dispersion a of 
the log-normal function is related to average gas density of the disk. We then calculate SFR as a function of critical 
density of loca l star formation and star star formation efficiency. This is a generalized version of the Schmidt-law 
(|Schmidtlll959l) . and it can be applied to various situations. 

An alternative way to study the global SFR in galaxies is simulating star formatio n directly in numerical models 
(e.g. iLi. Mac Low, fe Klessenl I^OllKravtsov fe Gnedinl [20051 : ITasker fe Brvar] 12001) . In this approach, 'stars' are 
formed according to a 'star formation recipe', and the resultant SFR is compared to observations. We do not take 
this methodology, because numerical modeling star formation in simulations still requires many free-parameters and 
assumptions. Moreover if the numerically obtained SFR deviates from the observed scaling relation (this is usually 
the case. See for example Tasker & Bryan 2006), it is hard to say what we can learn from the results. This deviation 
might be due to wrong implementation of star formation in the numerical code, or estimate of SFR in observations 
might be wrong, since SFR is not directly observable. One should also realize that comparison with observations of 
local galaxies is not necessarily useful when we discuss SFR in different situations, such as galaxy formation. In this 
paper, we avoid the ambiguity in terms of star formation in simulations, and alternatively discuss SFR based on an 
intrinsic statistical feature of the ISM. Effect of energy feedback from supernovae is complementarily discussed. 

This paper is organized as follows. In §2, we describe results of numerical simulations of the ISM in a galactic 
disk, focusing on the PDF. In §3, we start from a 'working-hypothesis', that is the inhomogeneous ISM, which is 
formed through non-linear development of density fluctuation, is characterized by a log-normal density PDF. After 
summarizing basic properties of the LN-PDF, we use the LN-PDF to estimate a global star formation rate, and it is 
compared with observations. In §4, we discuss implications of the results in §2 and §3. In order to distinguish local 
(i.e. sub-pc scale) and global (i.e. galactic scale) phenomena, we basically use cm -3 for local number density and M & 
pc~ 3 (or Mq pc -2 ) for global density 4 . 

2. GLOBAL SIMULATIONS OF THE ISM IN GALACTIC DISKS AND THE PDF 
2.1. Numerical Methods 

Evolution of rotationally supported gas disks in a fixed (i.e. time-independent), spherical galactic potential is 
investigated using three-dimensional hydrodynamic simulations. We take into account self-gravity of the gas, radiative 
cooling and heating processes. The nume rical scheme is an Euler method with a uniform Cartesian grid, which is 
based-on the code described in WN01 and IWadal (|2001| ). Here we briefly summarize them. We solve the following 
conservation equations and Poisson equation in three dimensions. 

^ + v-(H=o, (i) 

— + (-!;■ V)-u + -^ + V$ ox t + V$ sg = 0, (2) 

^2 + V • [(pE + p)v] = P T VV + r SN - p 2 A(T ff ), (3) 
V 2 $ sg = 4^G j0 , (4) 

where, p,p, v are the density, pressure, velocity of the gas. The specific total energy E = \v\ 2 /2+p/ (7 — l)p, with 7 = 
5/3. The spherical potential is $ext = -(27/4) 1 / 2 [« 2 /(r 2 +a 2 ) 1 / 2 +t- 2 /(r 2 +a 2 ) 1 / 2 1, where a 1 = 0.3 kpc a 2 = 5 kpc, and 
vi = v 2 = 200 km s" 1 . We also assume a cooling function A(T g ) (10 < T g < 10 8 K) (jSpaans fc Normarill997h with solar 
metallicity. We assume photo-electric heating by dust and a uniform UV radiation field, Tuv = l.Ox 10 _23 eG?o ergs s -1 , 
where the heating efficiency e is assumed to be 0.05 and Go is the incident FUV field normalized to the local interstellar 
value (jGerritsen fc Ickd 119971 ). In order to focus on a intrinsic inhomogeneity in the ISM due to gravitational and 
thermal instability, first we do not include energy feedback from supernovae. In §4.5 we show a model with energy 
feedback from supernovae (i.e. TgN 7^ 0). Note that even if there is no ra ndom energy input from supernov ae, turbulent 
motion can be maintained in the multi-phase, inhomogeneous gas disk f Wada. Meurer. fc No rman 2003). 

The hydrodynamic part of the basic equations is solved by AUSM (Advection Upstream Splitting Method) 
(jLiou fc St cffcn 1993) with a uniform Cartesian grid. We achieve third -order spatial accuracy with MUSCL (Monotone 
Upstream-centered Schemes for Conservation Laws) (jvan Le^[l97l . We use 512 x 512 x 64 grid points covering a 

1q pc — 3 ~ 66.7 cm -3 , if the mean weight of a particle is 0.61mjj- 



Wada, Norman 



3 



TABLE 1 

Initial density and fitting parameters of the density PDF 



Model 


ft a (Mq pc" 3 ) 


<no b 


a c 


log p d 


a c 


a"10,M* 


CFM g 




A 


5 


1.025 


2.360 


-1.40 


0.09 


1.025 


2.360 


2.483 


B 


10 


1.188 


2.735 


-1.55 


0.12 


1.188 


2.735 


2.769 


C 


15 


1.223 


2.816 


-1.60 


0.18 


1.273 


2.849 


2.810 


D 


50 


1.308 


3.012 


-1.50 


0.20 


1.308 


3.012 


3.104 



2.56 x 2.56 x 0.32 kpc 3 region (i.e. the spatial resolution is 5 pc). For comparison, we also run models with a 10 pc 
resolution. The Poisson equa tion is solved to calculate se lf-gravity of the gas using the Fast Fourier Transform (FFT) 
and the convolution method (jHocknev fc EastwoodllT98ll ). In order to calculate the isolated gravitational potential of 
the gas, FFT is performed for a working region of 1024 x 1024 x 128 grid points (see details in WN01). We adopt 
implicit time integration for the cooling term. We set the minimum temperature, for which the Jeans instability can 
be resolved with a grid size A, namely T m ; n = 35 K (p/10M Q pc _3 )(A/5 pc) 2 . 

The initial condition is an axisymmetric and rotationally supported thin disk (scale height is 10 pc) with a uniform 
density pi. We run four models with different pi (see Table 1). Random density and temperature fluctuations are 
added to the initial disk. These fluctuations are less than 1 % of the unperturbed values and have an approximately 
white noise distribution. The initial temperature is set to 10 4 K over the whole region. At the boundaries, all physical 
quantities remain at their initial values during the calculations. 

2.2. Numerical Results 

Figure [T] is snapshots of density distribution at a quasi-stable state for models A, B and D on the x-y plane and x-z 
plane. Depending on the initial density (pi = 5, 10, and 50 Mq pc -3 ), distribution of the gas in a quasi-equilibrium is 
different. In the most massive disk (Fig. lc, model D), the disk is fragmented into clumps and filaments, on the other 
hand, the less massive disk (Fig. la, model A) shows more axi-symmetric distribution with tightly winding spirals 
and filaments. In Fig. [2j we show vertical structures of density of model D (t = 42 Myr). The high density disk 
(p > 0.1Af Q pc~ 3 ) is resolved by about 10 grid points out of 64 total grid points for z-direction. Density change in 
about 5 orders of magnitude from the disk to halo gas is resolved. 

Figure [3] is time evolution of a density PDF in model D. The initial uniform density distribution, which forms the 
peak around p ~ 50M© pc -3 , is smoothed out in ~ 10 Myr, and it turns out to a smooth distribution in ~ 20 Myr. 
Figure 2] is PDFs spatially de-convolved to four components (thin disk, thick disk in the inner disk, halo, and the 
whole computational box. See figure caption for definitions) of model D in a quasi-steady state (t = 32 Myr). It is 
clear that dense regions (p > 0.01M© pc~ 3 ) can be fitted by a single log-normal function, /ln, over nearly 6 decades: 



/ln(/?; Po, o-io) = 



27rerio ln 10 



exp 



log'O/Po) 2 



2a 2 



(5) 



where a is volume fraction of the high-density part which is fitted by the LN-PDF. On a galactic scale, the ISM is 
a multi-phase, and dense, cold gases occupy smaller volumes than diffuse, hot gases. Therefore, it is reasonable that 
the density PDF shows a negative slope against density as shown in the previous numerical simulations. However, 
our results suggest that structure of t he ISM is not scale-free. The LN-PDF implies that formation process for the 
high density part is highly non- linear (|Vazquez-Semadenilll994[) . High density regions can be formed by convergent 
processes, such as mergers or collisions between clumps/filaments, compressions by sound waves or shock waves. Tidal 
interaction between clumps and galactic shear or local turbulent motion can also change their density structure. For 
a large enough volume, and for a long enough period, these processes can be regarded as many random, independent 
processes in a galactic disk. In this situation, the density in a small volume is determined by a large number of 
independent random events, which can be expressed by p = Itl^i 8ip s with independent random factors 8i and initial 
density p s . Therefore the distribution function of log(p) should be Gaussian according to the central limit theorem, if 
N -> oo. 

In Table 1, we summarize fitting results for the PDF in the four models in terms of the whole volume. One of 
important results is that there is a clear trend that the dispersion is larger for more massive systems. The less massive 
model (model B) also shows a LN-PDF (Fig. [5J, but the dispersion is smaller (ctio = 1.188) than that in model D. 
The PDF for the whole volume is log-normal, but there is a peak around logp = 1.3. As shown by the dotted line, 
this peak comes from the gas just above the disk plane (10 < z < 20 pc) in the inner disk (r < 0.38 kpc), where the 
gas is dynamically stable, and therefore the density is not very different from its initial value (pi — 10). This peak is 
not seen in the dotted line in model D (Fig. [J), since the density field in model D is not uniform even in the inner 
region due to the high initial density. 

The excess of the volume at low density (< 10 _3 Af Q pc -3 , see Figs. [4] and [5]) is due to a smooth component that 
is extened vertically outside the thin, dense disk. About 10-20% of the whole computational box is in the LN regime 
(i.e. a = 0.1 — 0.2). However, most fraction of mass is in the LN regime (see Fig. [3]). 

In order to ensure that gravitational instability in high density gas is resolved, we set the minimum temperature 
7min, depending on a grid-size (§2.1). We run model C with two different resolutions, 5 pc and 10 pc. As seen in 
Fig. [51 although the "tangled web" structure is qualitatively similar between the two models, typical scales of the 
inhomogeneity, e.g. width of the filaments, are different. The scale height of the disk H in the model with A = 5 pc 
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is about a factor of two smaller than the model with A = 10 pc. This is reasonable, since T min oc A 2 and H oc T 3 1/2 ,if 
the disk is vertically in a hydrostatic equilibrium. In Fig. we compare density PDF in the two models with different 
spatial resolution. PDFs in the models are qualitatively similar, in a sense that the PDF has a LN part for high density 
gas, although the PDF in the model with a 10 pc resolution is not perfectly fitted by a LN function. This comparison 
implies that the spatial resolution should be at least 5 pc in order to discuss PDF of the ISM in galactic disks. 

Figure [5] shows a phase-diagram of model D. Two dominant phases in volume, i.e. warm gas T g ~ 8000 K and cold 
gas (T g ~ 30 — 1000K) exist. Temperature of the gas denser than p ~ 10Af© pc~ 3 is limitted by T min . As suggested 
by Fig. [8j frequency distribution of temperature is not represented by a single smooth function, which is a notable 
difference to density PDF. High density gas (p > 10~ 2 M@ pc -3 ) is not isothermal, suggesting that isothermality is 
not a necessary condition for the LN density PDF in global models of the ISM (see also §4.1). 

As qualitatively seen in Fig. [1] the dispersion of LN-PDF is larger in more massive system. This is more quantitatively 
seen in Fig. [9l mass- weighted PDFs in 4 models with different pi. In the fit, we use the same a obtained from the 
volume-weighted PDF, and the mass-weighted characteristic density po,Af is calculated using 



Po.m = Poe" ■ (6) 

Equation ([6]) is always true if the PDF is log-normal. The dispersion a is 2.36 in model A (pi = 5M© pc~ 3 ) and 3.01 
in model D (pi = 50Af Q pc -3 ) (Table 1). The dependence of a on the initial gas density (or total gas mass) is a 
natural consequence of LN-PDF (see §3.1). From comparison between the volume-weighted and mass-weighted PDF, 
it is clear that even if the smooth non-LN regime occupy the most part of the volume of the galactic disk-halo region 
(i.e. a < 1), the mass of ISM is dominated by the LN-PDF regime. 

Another important point in the numerical results is that the characteristic density of LN-PDF, po, does not signif- 
icantly change among the models, despite the wide (almost 8 orders of magnitude) density range in the quasi-steady 
state. In fact, as shown in Table 1, po is in a range of 1.7 < po < 2.7 cm -3 among the models with pi = 5 — 50M@ 
pc" 3 . Physical origin of this feature is discussed in §3.1. 

In summary, using our three-dimensional high resolution hydrodynamic simulations, we find that the ISM in a 
galactic disk is inhomogeneous on a local scale, and it is consisted of many filaments, clumps, and low-density voids, 
which are in a quasi-steady state on a global scale. Statistical structure of the density field in the galactic disks is well 
described by a single LN-PDF over 6 decades in a high density part (i.e. p > O.OIM© pc -3 ). Most part of the mass is 
in the regime of LN-PDF. The dispersion of the LN-PDF is larger for more massive system. 

3. PROBABILITY DISTRIBUTION FUNCTION AND THE STAR FORMATION RATE 

3.1. Basic Properties of the Log-Normal PDF 

In this section, based on properties of LN-PDF, we discuss how we can understand the numerical results in §2. 
Suppose that the density PDF, f(p) in the galactic disk is described by a single log-normal function: 

Hp/po) 2 



f(p)d P 



1 



27TCT 



exp 



2a 2 



d\np, 



(7) 



where po is the characteristic density and a is the dispersion. The volume average density (p)v for the gas described 
by a the LN-PDF is then 



(Ph 



1 



27TO- 

= P.^I\ 





[ Hp/po) 2 ] 


/ pexp 


2d 2 



d In p, 



The mass-average density is 



(P) 



M - 



Jp 2 dV 



Poe 



where M t = J pdV is the total mass in the log-normal regime. The dispersion a of the LN-PDF is therefore 



In 



(p) 



M 



Equivalently, using p , 



:21n 



hi 



(P)v 

f {p)v 

V Po 
( (p)m 

V Po 



(8) 
(9) 

(10) 

(11) 

(12) 
(13) 



Suppose the characteristic density po is nearly constant as suggested by the numerical results, Equations (fT2")) or (|13p 
tells that the dispersion a is larger for more massive systems, which is also consistent with the numerical results. For 
a stable, uniform system, i.e. (p)y = Pen o~ should be zero. In another extreme case, namely (p)v ~> oo, ~ > °o, but 
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this is not the case, because the system itself is no longer dynamically stable. Therefore, a should take a number in an 
appropriate range in a globally stable, inhomogeneous systems. Galactic disks are typical examples of such systems. 

In our numerical results, the ISM with low density part (typically less than 10 _3 M© pc~ 3 ) is not fitted by the 
LN-PDF. If density field of the ISM in a fraction, a, of the arbitrarily volume is characterized by the LN-PDF, the 
volume average density in the volume is p — a{p)v- In this case, eq. (|12[) is modified to 

o* = 2)n(-M. (14) 

As shown in Table 1, a in the numerical results well agrees with a p in each model. In numerical simulations, it is 
easy to know the volume in the LN-regime (i.e. ax volume of the computational box), therefore calculating a p is 
straightforward using eq. (|14[) . However, if one wants to evaluate a from observations of galaxies, it is more practical 
to use the mass-average density, cq. ([13")) , because it is expected that mass of the galactic ISM is dominated by the 
LN-regime in mass (see Fig. 

Interestingly, even if the density contrast is extremely large (e.g. 10 6 — 10 8 ), the numerical results show that the 
characteristic density po is not sensitive for the total gas mass in a kpc galactic disk (po = 1-7 — 2.7 cm -3 ). If this 
is the case, what determines po? In Fig. [101 we plot effective pressure, p e s = Pth + pv 2 as a function of density, 
where p t h is the thermal pressure, and v t is turbulent velocity dispersion, which is obtained by averaging the velocity 
field in each sub-region with a (10 pc) 3 volume. If the turbulent motion in a volume with a size L is originated in 
selfgravity of the gas and galactic rotation iWada. Meurer. fc Normanl [2002:) , p c g ~ pGM g /L cx p 4 / 3 , provided that 
the mass in the volume M g is conserved (i.e. p cx Zr 3 ) 1 . In fact, Fig. [TOl shows that a majority of the gas in the 
dense part (p > 10 _1 M Q pc -3 ) follows p cr f cx p 4 / 3 . On the other hand, for the low density regime, p c s cx p 2 / 3 . In the 
present model, a dominant heating source in low density gases is shock heating. Shocks are ubiquitously generated by 
turbulent motion, whose velocity is ~ several 10 km s . The kinetic energy is thermalized at shocks, then temperature 
of low density gas goes to T g ~ 10 4-5 K. In supersonic, compressible turbulence, its energy spectrum is expected to be 
E{k) cx fc~ 2 , where total kinetic energy E t = J E(k)dk. Therefore p c g cx pv 2 cx p 2 / 3 , because v 2 cx fc -1 cx p -1 / 3 . Then 
the 'effective' sound velocity c 2 ff = dp e e/dp cx p" 1 ! 3 and c 2 ff cx p 1 ! 3 for the low and high density regions, respectively. 
This means that the effective sound velocity c e ff has a minimum at the transition density (pt) between the two regimes. 
In other words, since c 2 ff = c 2 + v 2 , there is a characteristic density below/above which thermal pressure/turbulent 
pressure dominates the total pressure, i.e. c 2 ~ v 2 . Thus, 

^ ~ ^ ~ Gp t L 2 , (15) 

p Lt 

where p is average mass per particle, and L is size of the largest eddy of gravity-driven turbulence, which is roughly 
the scale height of the disk. As seen in Fig. gas temperature around po = 10 -1 ' 5 Mq pc" 3 is T g ~ 80 K. The scale 
height of the dense disk is about 10 pc (Fig. [2]) in the present model, therefore 

This transition density is close to the characteristic density in the numerical results, i.e. po = 10~ 14 — 10 _16 M© 
pc~ 3 (« 2.7 - 1.7 cm -3 ). In a high-density region (p > p t ~ po), stochastic nature of the system, which is the origin of 
the log-normality, is caused in the gravity- driven turbulence, and for low density regions its density field is randomized 
by thermal motion in hot gases. For the gases around po, the random motion is relatively static, therefore probability 
of the density change is small. As a result, the density PDF takes the maximum around po- An analogy of this 
phenomenon is a snow or sand drift in a turbulent air. The material in a turbulent flow is stagnated in a relatively 
"static" region. 

As seen in Fig. [2j the present gas disk has a much smaller scale height than the ISM in real galaxies, which is 
about 100 pc. If turbulent motion in galactic disks is mainly caused by gravitational and thermal instabilities in a 
rotating disk (jWada. Meurer. fc Norman]|2002f) . we can estimate p t by eq. (Ti"rJ)) in galactic di sks, but gas temperatur e 
of the dominant phase in volume is T g ~ 10 4 K, as suggested by two-dimensional models (W ada Sz Normanl l2001h . 
Therefore, we expect that pt ~ 1 cm~ 3 is also the case in real galaxies. However, if the turbulent motion in ISM is 
caused by different mechanisms, such as supernova explosions, magneto-rotational instabilities, etc., dependence of 
Pes on density could be different from p e g cx p 4 ^ 3 , and as a result pt ~ const. ~ 0(1) cm -3 could not be always true. 
Determining po in a much larger disk than the present model, taking into account various physics is an important 
problem for 3-D simulations with a large dynamic range in the near future. In §3.3, we compare our results with SFR 
in spiral and starburst galaxies, in which po is treated as one of free parameters (see Fig. [12] and related discussion). 

Suppose the volume average density is (p)y = 3A/ Q pc~ 3 , we can estimate the dispersion a is a ~ 3 for p = 1 cm~ 3 
from eq. (JT2J) . For a less massive system, e.g. (p)v = O.3M pc~ 3 , a ~ 2.1. Therefore we can expect a larger star 
formation rate for more massive system, since a fraction of high density gas is larger. Using this dependence of a on 
the average gas density, we evaluate global star formation rate in the next section. 

1 L is not a size of molecular clouds, but size of 'eddies' turbulent motion in the inhomogeneous ISM. (See Fig. 2 in 
Wada, Meurer. & Norman 2002, for example). Therefore, the Larson's law, i.e. p oc (size of molecular clouds) -1 , is not a relevant 
scaling relation here. 
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3.2. Global Star Formation Rate 

Based on the numerical results in §2 and the properties of LN-PDF described in §3.1, we here propose a simple 
theoretical model of the star formation on a glob al scale. We ass ume that the multi-phase, inhomogeneous ISM in a 
galactic disk, can be represented by a LN-PDF. [Elmcgrccn] (|2002f) discussed a fraction of high density gas and global 
star formation rate assuming a LN-PDF. Here we follow his argument more specifically. 

If star formation is led by gravitational collapse of high-density clumps with density p Cl the star formation rate per 
unit volume, /?'*, can be written as 

p\ = e c {G Pc ) l / 2 f c {p) v (17) 

where / c (p c , er) is a mass fraction of the gas whose density is higher than a critical density for star formation (p > p c ), 
and e c is the efficiency of the star formation. 

Suppose the ISM model found in §2, whose density field is characterized by LN-PDF, is applicable to the ISM in 
galactic disk, f c is 



fc(S c , er) = 





[ 2a 3 


dQnS) 




(ln<5) 2 " 
2cr 2 


d(ln«S) 



\(l-Ed[z(S c ,a)]) 



where 8 = p/p , S c = p c /po, and 



z(S c ,a) 



In 5r 



V2a 



(18) 



(19) 



(20) 



The fraction of dense gas, f c is a monotonic function of 5 C and a, and it decreases rapidly for decreasing er. Suppose 
S c = 10 5 , f c ~ 10~ 2 for a = 3.0, and f c ~ 10~ 6 for a = 2.0 2 . The star formation rate (eq. (JTTj) ) per unit volume then 
can be rewrite as a function of e c , 5 C , and er: 



p i ,(e c ,S c ,a) = e c (G6 c ) 1 / 2 f c p s /2 e^ 2 
Using equations (TT2")) . (fT!))) and (j2"0)) , we can write SFR as a function of volume- average density (p)v> 



(Ph 



lM©pc 



-3 



Pa 



1 cm 



-3 



10 5 < 



3.6 x 10 _7 e c M yr" 



^pc- 3 



(21) 



(22) 



1 - Erf 



HPcPo/(p)y) 
2[ln«p) y /po)] 1 / 2 



Pl /2 (p)v- 



We plot eq. ([22l as a function of the volume average density in Fig. [TTJ Four curves are plotted for p c = 10 3 , 10 4 , 10 5 
and 10 6 cm -3 . We can learn several features from this plot. SFR increases rapidly for increasing the average density, 
especially for lower average density (< 100 A/© pc" 3 ) and higher critical density. This behavior is naturally expected 
for star formation in the ISM described by LN-PDF, because the dispersion a changes logarithmically for the average 
gas density (eq. |12|)). For large density, it approaches to SFR oc (p)y SFR does not strongly depend on p c around 



(Ph 
to p, 



1M Q pc , and SFR is larger for larger p c beyond the density. This is because the free fall time is proportional 



-1/2 



and this compensates decreasing the fraction of high density gas, f c 



Fig. [T2l shows how SFR changes as a function of po- SFR does not strongly depend on p around P q < O.IMq pc 3 
(~ 7 cm" 3 ). See also Fig. rjj]in terms of comparison with observed SFR. 

3.3. Comparison with Observations 

The theoretical SFR based on the LN-PDF in the previous section has a couple of free parameters. Comparison 
between the model and observations is useful to narrow down the parameter ranges. 

We should note, however that the gas disk presented in §2 is geometrically thiner (~ 10 pc) than the real galactic 
disk (~ 100 pc). In this sense, although the present models are based on full 3-D simulations, they are not necessarily 
adequate for modeling typical spiral galaxies 3 . There are couple of reasons why the present model has small scale 
height. Since the disks here are relatively small (r ~ 1 kpc), the disks tend to be thin due to the deep gravitational 
potential. Self-gravity of the gas also contributes to make the disks thin. Radiative energy loss in the high density 
gas in the central region cancels the energy feedback from the supernovae (see §4.5). Therefore, a simple way to make 
the model gas disks thicker to fit real galactic disks is simulating a larger disk (e.g. r ~ 10 kpc) with an appropriate 



2 In I Elmcgrccn (2002), the dispersion of the LN-PDF, a = 2.4 is assumed, which is taken from 2-D hydrodynamic simulations of the 
multi-phase ISM in WN01. 

3 One can refer to a critical study for modeling galactic ISM using a 2-D approximation ( Sanchcz-Salccdo 2001) 
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galactic potential and supernova feedback, solving the same basic equations 4 . Besides supernova explosions, there 
are possible physical processes to puff up the disks, such as nonlinear development of magneto-rotational instability, 
heating due to stellar wind and strong radiation field from star forming regions. These effects should be taken into 
account in 3-D simulations with high resolutions and large dynamic ranges in the near future, and it is an important 
subject whether the LN-PDF is reproduced in such more realistic situations. Another interesting issue in terms of PDF 
is effect of a galactic spiral potential. However, here we try to make a comparison with observations, assuming that 
the log-normal nature in the density field found in our simulations is also the case in real galactic disks. This would 
not be unreasonable, because log-nor mality is independent of th e geometry of the system. In fact a LN-PDF is also 
found in a 3-D torus model for AGN (jWada. fc To misaka 2003). A LN-PDF is naturally expected, if the additional 
physical processes causes non-linear, random, and independent events that change the density field. 

Figure [T3l shows su rface star formation rate E* (M© yr" 1 pc~ 2 ) in normal and starburst galaxies (Komugi ct al. 
I2005t iKennicuttl fl998l ) . The scale height of the ISM is assumed to be 100 pc. Four curves represent model SFR with 
different critical density. It is clear that smaller 8 C is preferable to explain the observations, especially for low average 
density. For example, SFR with 8 C = 10 5 is too steep. Similarly, Fig. [l4]shows dependence of SFR on the characteristic 
density pq. As mentioned above, SFR is not very sensitive for changing po, but po < 1.0 cm -3 is better to explain the 
observations. From Fig. 1151 which is how SFR depends on a fraction of LN part in volume, a. SFR is sensitive for 
a especially for low density media. The plot implies that a = 0.01 is too steep to fit the observations, especially for 
normal galaxies. 

After exploring these parameters, we find that Fig. [16] is the best-fit model with po = 1 cm~ 3 , a = 0.1 (eq. fl3])) 
and 8 C = 10 3 . Four curves corresponds to SFR with efficiency, e c = 0.1, 0.01, and 0.001. The model slope approaches 
to SFR cx Ti g for large density, which is shallower than the Kennicutt law (i.e. SFR cx S^ 41 ). Observed SFRs in most 
galaxies are distributed between the model curves with e e ~ 0.1 and 0.001 . The starburst galaxies are distributed in 
e c = 0.1 — 0.01, on the other hand, the normal galaxies (jKomugi et al.ll2"005h show systematically smaller SFR, which 
is consistent with smaller efficiency (i.e e c = 0.01 — 0.001) than those in starburst galaxies. This suggests that the large 
SFR in starburst galaxies is achieved by both high average gas density and large (several %) star formation efficiency 
i n dense clouds 5 . 

iGao fc Solomon! (|2004f ) found that there is a clear positive correlation between HCN and CO luminosity in 65 infrared 
luminous galaxies and normal spiral galaxies. They also show that luminous and ultra-luminous infrared galaxies tend 
to show more HCN luminous, suggesting a larger fraction of dense molecular gas in active star forming galaxies. Figure 
[171 is SFR as a function of f c using the model in the previ ous section. Interestingly, a qualitative trend of this plot 
is quite similar to that of Fig. 4 in IGao fc Solomon] (|2004f ) . which is SFR as a function of Lhcn/£co(°< dense gas 
fraction). Both observations and our model show that SFR increases very rapidly for increasing f c , especially when 
f c or Lhcn/Lco is small (< 0.1) and then it increases with a power-law. 

4. DISCUSSION 
4.1. What Determines the Dispersion of LN-PDF? 

Elmegreen (2002) first pointed out that if the density PDF in the ISM is described by a log-normal function, the star 
formation rate should be a function of the critical density for local star formation and the dispersion of the LN-PDF. 
The critical density for local star formation should be determined not only by gravitational and thermal instabilities of 
the ISM, but also by magneto- hydrodynamical, chemical and radiative processes on a pc/sub-pc scale. It is beyond a 
scope of the present paper that how the critical density is determined. The other important parameter, a in LN-PDF, 
should be related to physics on a global scale. 

Some authors claimed that the LN-PDF is a char acteristic feature in an i s othermal, turbulent flow, and its 
dispersion is determined by the rms Ma ch number (|Vazquez-Semadenl 11994 iPadoan. Nordlund. fc JohnU Il997t 
iNordlund fc Padoanlfl998tlScalo et al.lll~998h . This argument might be correct for the ISM on a lo cal scale, for example 
an int ernal structure of a giant molecular cloud, which is nearly isothermal and a single phase. iKrumholz fc McKed 
(2005) derived an analytic prediction for the star formation rate assuming that star formation occurs in virialized 
molecular clouds that are supersonically turbulent and that density distribution within a cloud is log-normal. In this 
sense, their study is similar to Elmegreen (2002) and the present work. However, they assumed that the dispersion of 

the LN-PDF is a function of "one-dimensional Mach number", A4, of the turbulent motion, i.e. e CT w 1 + 3A^ 2 /4, 
which is suggested by numerical experiments of isothermal turbulence (Padoan & Nordlund 2002). A similar empirical 
relation between the density contrast and magneto-sonic Mach number is suggested by Ostriker, Stone & Gammie 
(2001). Yet the physical reason why the dispersion depends on the Mach number is not clear. If high density regions 
are formed mainly through shock compression in a system with the rms Mach number M rms , the average density 
contrast (Sp/p) ~ (p)v/po — 1 could be described by M Tms , and it is expected that (p)v/po — e a / 2 = 1 + A4 2 ms using 
eq. ([9]). However, the ISM is not isothermal on a global scale. An inhomogeneous galactic disk is characterized by 

4 This is practically difficult, if a pc-scale spatial resolution is required. A recent study by Taskcr & Bryan (2006) is almost the only 
work that can be directly compared with real galactic disks. Unfortunately, Taskcr & Bryan (2006f) do not discuss on the density PDF, 
but they take a different approach to study SFR by generating "star" particles (see §1). They claimed that a Schmidt-law type SFR is 
reproduced in their simulations. This is consistent with our analysis, provided that the density PDF is log-normal like, and that local star 
formation takes place above a critical density. 

° An alternative explanation can be possible: the typical scale height of the star forming regions is different by a factor of 10-100 in 
normal and starburst galaxies for the same efficiency. 
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a fully developed turbulence, and its velocity dispersion is a function of scales as shown by power-law energy spectra 
(jWada. Meurer. fe Normanll2002l ) . This means that the galactic disk cannot be modeled with a single 'sound velocity' 
or velocity dispersion of the gas. Therefore we cannot use the empirical relation on a for galactic disks. 

One should note that there is another problem on the argument based on the rms-Mach number, which is in terms 
of origin and mechanism of maintaining the turbulence in molecular cl ouds. Numerica l experiments suggest ed that the 
turbulence in molecular clouds is dissipated in a sound crossing time (|Mac Lowlll999t [Ostriker et al.ll2001l ). and there 
is no confirmed prescriptions on energy sources to compensate the dissipation. Therefore, it is more natural to assume 
that the velocity dispersion in a molecular cloud is not constant in a galactic rotational period, and as a result the 
structure of the density field is no longer static. If this is the case, taking the rms-Mach number as a major parameter 
to describe the global star formation rate would not be adequate. The decaying turbulence could be supported by 
energy input by supernovae or out-flows by proto-stars, but even in that case there is no clear reason to assume a 
uniform and time-independent Mach number. 

4.2. Observational evidence for LN-PDF 

It is practically difficult to know the PDF of the ISM in galactic disks directly from observations. We have to map 
the ISM in external galaxies with fine enough spatial resolution by observational probes that cover a wide density 
range. The Atacama Large Millimetcr/Submillimcter Array will be an ideal tool for such observations. Nevertheless, 
there is indirect ev idence to support the LN-PDF in the Large Magellanic Cloud ( LMC), which is mapped by HI with 
a 15 pc resolution (Kim et al. 2003) and by CO ( J = 1 — 0) with 8 pc resolution (|Fukui et al.ll200it lYamaguchi et al.l 
2001). We found that the distribution function of HI intensity distribution and a mass spectrum of CO clouds are 
consistent with a nu merical model of the LMC, in which the density PDF in the simulation is nicely fitted by a single 
log-normal function (|Wada. Spaans, fc Kiml 12000). Although we need more information about density field by other 
probes, this suggests that the entire density filed of the ISM in LMC could be modeled by a LN-PDF. 

Recently, iTozzi. et al.1 ()2006f ) claimed that distribution of absorption column density N-r in 82 X-ray bright sources 
in the Chandra Deep Field South is well fitted by a log-normal function. This suggests that the obscuring ma- 
terial around the AGNs is inhomogeneous as suggested by previous numerical simulations (|Wada fc Normanl 120021 : 
IWada. fc~T omisaka] |2005[ ) . and their density field is log-normal, if orientation of the obscuring "torus" is randomly 
distributed for the line of sight in the samples. 

4.3. The Q criterion for star formation 

The origin of the Schmidt-Kennicutt relation on the star formation rate in galaxies (|Kennicuttl Il998f ) has been 
often discussed in terms of gravitational instability in galactic disks. More specifically it is claimed that the thresh- 
old density for star formation can be represented by the density for which the Toomre Q parameter is unity, i.e. 
Q = kc s /(ttGY, ) = 1, where k, c s , and E g are epicyclic frequency, sound velocity, and surface density of gas (|Kennicuttl 
ll998HMartin fc Kennicuttll2001h. However , this is not supported by recent observations in some gas-rich spiral galax- 
ies (jWong fc Blitz! 120021 : iKoda et al.l [2005h . Based on our picture described in this paper, it is not surprising that 
observationally determined Q or the critical density do not correlate with the observed star formation rate. Stars 
are formed in dense molecular clouds, which is gravitationally unstable on a local scale (e.g. < 1 pc), but this is 
basically independent of the global stability of the galactic disk. Even if a galactic gas di sk is globally stable (e.g. 
effective Q > 1), cold, dense molecular clouds should exist (see also lWada fc Normanl 1 19991 ). Numerical results show 
that once the galactic disk is gravitationally and thermally unstable, inhomogeneous structures are developed, and 
in a non-linear phase, it becomes "globally stable" , in which the ISM is turbulent and multi-phase, and its density 
field is characterized by LN-PDF. In that regime, SFR is determined by a fraction of dense clouds and by a critical 
density for "local star formation", which should be related to physical/chemical conditions in molecular clouds, not 
in the galactic disk. In this picture, SFR naturally drops for less massive system without introducing critical density 
(see Figs. E] and [HI). 

Another important point, but it has been often ignored, is that the Q criterion is de rived from a dispersion re lation 
for a tightly wound spiral perturbation in a thin disk with a uniform density (see e.g. iBinnev fc Tremainl ll987). The 
equation of state is simply assumed as p oc S 7 , where £ is surface density and 7 is a constant. The criterion for 
instability, i.e. Q = 1 means that the disk is linearly unstable for an axisymmetric perturbation. These assumptions 
are far from describing star formation criterion in an inhomogeneous, multi-phase ga lactic disks. One should also no te 
that gas disks could be unstable for non-axisymmetric perturbations, even if Q > 1 dGoldreich fc Lvnden-Belll[l965h . 

Finally, we should emphasize that it is observationally difficult to determine Q precisely. All the three variables in 
the definition of Q, i.e. surface density of the gas, epicyclic frequency, and sound velocity, are not intrinsically free from 
large observational errors. Especially, it is not straightforward to define the 'sound velocity', c s in the multi-phase, 
inhomogeneous medium, and there is no reliable way to determine c s by observations. Therefore one should be careful 
for discussion based on the absolute value of Q in terms of star formation criterion. 

4.4. LN-PDF and star formation rate in simulations 

iKravtsov fc Gnedinl (p005h run a cosmological N-body/AMR simulation to study formation of globular clusters in 
a Milky Way-size galaxy. They found that the density PDF in the galaxy can be fitted by a log-normal function, and 
it evolves with the redshift. The dispersion for the exponent of a natural log is a ~ 1.3,2.0, and 2.8 at z = 7,3.3 
and 0, respectively. As shown in the previous section, the dispersion of the LN-PDF is related with the average gas 
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density in the system. Therefore, the increase of the dispersion in the galaxy formation simulation should be because 
of increase of the gas density and /or decrease of the characteristic density po. The former can be caused by accretion 
of the gas, in fact. IKravtsov fc Gned in (2005) shows that the gas density increases until z ~ 5. They also shows that 
the characteristic density po decreases from z ~ 7 to z ~ 3. Based on the argument on p e ff in §3.1, this is reasonable, if 
1) a clumpy, turbulent structure is developed due to gravitational instability and 2) temperature of the lower density 
gas decreases due to radiative cooling with time. Qualitatively this is expected because at lower redshift, heating 
due to star formation and shocks caused by mergers are less effe c tive, but radiative cooling becomes more efficient 
due to increase of metallicity and gas density. IKravtsov fc Gnedin! (|2005l ) mentioned that widening the LN-PDF with 
decreasing redshift is due to increase of the rms Mach number of gas clouds. However, we suspect that this is not the 
c ase in formation of galact i c disk s (see discussion in §4.1). 

iLi. Mac Low. &: K lcsscn ( 2006) investigated star formation in an isolated galaxy, using three-dimensional SPH code 
(GADGET). They use sink particles to directly measure the mass of gravitationally collapsing gas, a part of which 
is considered as newly formed stars. They claimed that the Schmidt law observed in disk galaxies is quantitatively 
reproduced. They suggest that the non-linear development of gravitational instability determines the local and global 
Schmidt laws and the star formation rate. Their model SFR shows a rapid decline for decreasing the surface density, 
which is consistent with our prediction (Fig. 11). 

4.5. Effects of energy feedback 

In §2, we focus on models without energy feedback from supernovae in order to know an intrinsic structures of 
the ISM, which is dominated by gravitational and thermal instabilities. However, one should note that the LN-PDF 
is robust for including energy feedback from supernovae as suggested by previous 2-D models for galactic disks and 
3-D models for galactic central regions (WN01 and W01). This is reasonable, because stochastic explosions in the 
inhomogeneous medium is a preferable situation for the LN-PDF. In order to confirm this in 3-D on a galactic scale, 
we run a model with energy feedback from supernovae, in which a relatively large supernova rate (1.5 x 10~ 3 yr _1 
kpc -2 ) is assumed. Energy from a supernova (10 51 ergs) is injected into one grid cell which is randomly selected in 
the disk. In Figs. fl8l and [T9l we show density distributions and PDFs in models with and without energy feedback. 
It is clear that even for the large supernova rate the density distribution and PDF are not significantly different from 
the model without energy feedback, especially for the regime above the characteristic density. 

4.6. Origin of a bias in galaxy formation 

It is observationally suggested that massive galaxies terminat e their star formation at higher redshift th an less 
massive galaxies, namely "down-sizing" in galaxy formation (e.g. iKauffma nrTet al.l [2003: Ko dama et al.|[2004f) . This 
seems to be inconsistent with a standard hierarchical clustering scenario. It is a puzzle why the star formation time 
scale is shorter in more massive galaxies, in other words, why the star formation rate and/or star formation efficiency 
are extremely biased in an environment to form massive galaxies. This might be understood by our result of the 
global star formation, that is SFR increases extensively as a function of average gas density (Fig. 11), e.g. SFR oc p 2 
(p ~ 1M pc -3 , for S c = 10 4 ). Since the average baryon density is higher in an environment where massive galaxies 
could be formed, SFR could be 100 times larger, if the average gas density is 10 times higher. It would be interesting 
to note that this tendency is more extreme for higher critical density (p c )- Critical density for s tar formation can 
be af fected by UV radiation through photoevapolation of dense clouds in proto-galactic halos (e.g. ISusa &: Umemural 
2004). Qualitatively this means that for a stronger UV field, the critical density is larger, therefore SFR can strongly 
depend on the average gas density. 

5. CONCLUSION 

Three-dimensional, high resolution hydrodynamic simulations of a galactic disk shows that the density probability 
distribution function (PDF) is well fitted by a single log-normal function over 6 decades. The dispersion of the log- 
normal PDF (LN-PDF), a can be described by a = \/2\n((p) m / Po) , where po is a characteristic density, and (p)m is 
mass- weighted average density of the ISM (eq. (I13|)). If all the ISM is in the regime of LN-PDF, the star formation rate 
can be represented by eq. . We find that star formation rate (SFR) is sensitive for increasing average gas density, 
especially for smaller (p)y and larger p c . It is however that SFR does not significantly depend on the characteristic 
density pq. We compare the observed star formation rate in normal and spiral galaxies, and find that a model with 
p c ~ 10 3 cm -3 , a volume fraction of the LN-part ~ 0.1, and po ~ 1 cm~ 3 , well explains the observed trend of SFR 
as a function of average surface density. If the scale height of the ISM in star forming regions is ~ 100 pc, the star 
formation efficiency in starburst galaxies is 0.1-0.01, and it is one-order smaller in normal galaxies. 

The log-normal nature of the density field should be intrinsic in an inhomogeneous, multi-phase ISM, if 1) the whole 
system is globally quasi-stable in a long enough period (for a galactic disk, it is at least a few rotational periods), 2) the 
system is consisted of many hierarchical sub-structures, and 3) density in such sub-structure is determined by random, 
non-linear, and independent processes. If these conditions are satisfied, any random and nonlinear processes that affect 
a density field should cause the log-normal probability distribution function. In this sense, most physical processes 
expected in real galactic disks, such as nonlinear development of magneto-rotational instability, interactions between 
the ISM and stellar wind, and heating due to non-uniform radiation fields originated in OB associations should also 
generate the log-normal PDF. These effects on the PDF could be verified in more realistic numerical simulations with 
a wide dynamic range in the near future. Observational verification on the density PDF of the ISM in various phases 
is also desirable. 
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Fig. 1. — Surface sections of the density distribution of the gas in models with different initial density (model A, B, and D). x-y and x-z 
planes are shown. Unit of length is kpc. (a) model A {pi = 5Mq pc — 3 , t = 103 Myr), (b) model B (pi = IOMq pc — 3 , t = 36 Myr), and 
(c) model D (p l = 50Af Q pc -3 , t = 32 Myr). 
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Fig. 2. — Vertical density structure of model D. Two lines represent density averaged for x-dircction (black dots) and y-direction (open 
circles). 




Log(density) 



Fig. 3. — Time evolution of density PDF for the whole computational volume for model D. The lines represent the PDF at t = 3, 10, 16, 21, 
and 43 Myr (thick solid line). The vertical axis is the number fraction of grid cells for the density. The unit of density is M© pc~ 3 . 
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Fig. 4. — Density PDF for four sub-regions in model D. Solid line is a thin disk (\z\ < 10 pc), dotted line is a inner disk (10 < z < 20 
pc and \x,y\ < 0.38 kpc), dashed line is a halo (30 < z < 160 pc), where z = is the galactic plane. The two solid curves are a log-normal 
function with <rio = 1.273, log(po) = —1.2 and a = 0.02 (disk), and <xio = 1.308, log(po) = —1.5 and a = 0.2 (the whole region). 
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Fig. 5. — Same as Fig. [1] but for model B. 
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Fig. 6. — Density structures of model C with two deferent resolutions. A temperature floor depending on the resolution: (a) A = 10 pc, 
and (b) A = 5 pc. T m ; n = 35 K (p^/IOMq pc — 3 ) (A/5 pc) 2 is assumed to ensure that the Jeans unstable is resolved. 
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Fig. 7. — Density PDF for model C with two deferent resolutions, (a) A = 10 pc, (b) A = 5 pc. Three lines are PDF for the whole 
computational box, disk, and the inner disk (dotted line). 
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Fig. 10. — Effective pressure-density diagram of model E. Two power-laws are p e ff oc p 2 / 3 and p c g oc p 4 / 3 . p c g = p tll + pi) 2 , where p t h 
is the thermal pressure and vt is turbulent velocity of the medium. The contours represent log-scaled volumes for given density and p e ff. 
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FlG. 11. — Star formation rate (p*) as a function of volume average density (p) for 0.01. po = 1 cm 3 is assumed. Four curves are 
plotted for <5 C = 10 3 (thick solid line), 10 4 , 10 s and 10 6 (from thick to thin lines) are plotted. 
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Fig. 12. — Star formation rate (p*) as a function of the characteristic density of the LN-PDF (po) for different average density 
(p)v = 10 2 , 10, and 1M Q pc~ 3 . <5 C = 10 4 is assumed. 
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Fig. 13. — Comparison between our models and observed surface star formation rate in ter ms of surface gas density. Black dots are 
starburst galaxies in Kennicutt (1998) and gray dots are normal galaxies in Komugi ct al. (2005). Four curves are SFR for S c = 10 2 (thick 
solid line), 10 3 , 10 4 and 10 5 . a = 0.1, po = 1-0, and e c = 0.01 are assumed. 
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Fig. 14. — Same as Fig. 1131 but for dependence of characteristic density of LN-PDF, po = 0.1 (thick solid line), 1 and IOMq pc" 
a = 0.1, <5 C = 10 3 , and e c = 0.01 are assumed. 
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Fig. 15. — Same as Fig. 1131 but for dependence of fraction of LN part, a = 0.01 (thick solid line), 0.1, and 1.0. po = 1.0, S c = 10 3 , and 
e c = 0.01 are assumed. 
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Fig. 18. — Same as Fig. 1, but comparison between models with and without energy feed back from supernovae. (a) model D, (b) model 
D, but with supernova explosions at a rate of 1.5xl0 — 3 yr — 1 kpc - 2 . 
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Fig. 19. — PDF in models with and without energy feedback from supernovae. Thin line is model D, and thick line is model D, but with 
supernova explosions at a rate of 1.5xl0 — 3 yr — 1 kpc -2 . 



